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ON THE ROBUSTNESS OF Z = 



1 GALAXY SIZE MEASUREMENTS THROUGH MODEL AND 
NON-PARAMETRIC FITS 
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Abstract 

We present the size-stellar mass relations of nearby {z = 0.01 — 0.02) SDSS galaxies, for samples 
selected by color, morphology, Sersic index n, and specific star formation rate. Several commonly- 
employed size measurement techniques are used, including single Sersic fits, two-component Sersic 
models and a non-parametric method. Through simple simulations we show that the non-parametric 
and two-component Sersic methods provide the most robust effective radius measurements, while 
those based on single Sersic profiles are often overestimates, especially for massive red/early- type 
galaxies. Using our robust sizes, we show that for all sub-samples, the mass-size relations are shallow 
at low stellar masses and steepen above ~ 3 — 4 x lO^^M©. The mass-size relations for galaxies 
classified as late-type, low-n, and star-forming are consistent with each other, while blue galaxies 
follow a somewhat steeper relation. The mass-size relations of early-type, high-n, red, and quiescent 
galaxies all agree with each other but are somewhat steeper at the high-mass end than previous 
results. To test potential systematics at high redshift, we artificially redshifted our sample (including 
surface brightness dimming and degraded resolution) to 2; = 1 and re-fit the galaxies using single 
Sersic profiles. The sizes of these galaxies before and after redshifting are consistent, and we conclude 
that systematic effects in sizes and the size-mass relation at z ~ 1 are negligible. Interestingly, since 
the poorer physical resolution at high redshift washes out bright galaxy substructures, single- Sersic 
fitting appears to provide more reliable and unbiased effective radius measurements at high z than 
for nearby, well-resolved galaxies. 
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1. INTRODUCTION 

Correlations among galaxy physical parameters such 
as stellar mass, luminosity, size, velocity dispersion, and 
their evolution with cosmic time are crucial for under- 
standing the formation and evolution of galaxies and im- 
posing constraints on theoretical models of their struc- 
tural assembly. Morphological scaling relations such as 
the relation between size and surface b rightness, t he cor- 
relation between size and luminosity (jKormendvl Il985) . 
and t he relation betwee n the effective radius and stellar 
mass (jShen et al.l I2003L hereafter SOS), vary for differ- 
ent types of galaxies. The differences between surface 
brightness profiles and sizes of galaxies are the products 
of the different physical processes governing their forma- 
tion and evolution. Precise measurements of these galaxy 
properties at low and high redshifts thus provide strong 
constraints on models of galaxy formation and evolution. 

Among these relations is the observed correlation be- 
tween half-light radius (size) and stellar mass, which 
is shown for the local Universe (S03) and persists 
up to very high r edshifts (e.g., IPaddi et aTl |2005|; 
Truiillo et ar '2006'; 'Fra nx et al.ri2 008>: ' Buitrago et"aLl 
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These authors also pointed 
out that sizes of galaxies at fixed stellar mass decrease 
as the redshift increases, i.e., galaxies were smaller in the 
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past. For instance, massive quiescent galaxies at z 2 
are about a facto r of ^ 6 sm aller tha n their c ounterp arts 
at z ^ (e.g., IPaddi et al. 2005; | van Do kkum et al.l 
I2OO8I ). Understanding the mechanism of the size evo- 
lution and how galaxies reach the mass-size relation at 
z = requires measuring these properties, especially 
sizes, very robustly in different redshift ranges. 

One of the main concerns is the accuracy of galaxy 
size determination at high redshifts. Galaxies at higher 
redshifts are at larger distances and therefore are dim- 
mer and have smaller apparent angular sizes. The low 
surface brightness envelopes of galaxies could fade away 
due to cosmological dimming and could have lower sig- 
nal to noise ratios (S/N), hence, may potentially bring 
systematics on their real size measurements. For exam- 
ple, the outer parts of early- type galaxies normally fade 
away gradually into the background sky noise, and it is 
very hard to define precise edges for these types of galax- 
ies. Underestimating the sizes of these galaxies at high 
redshifts co uld have an effect o n the inferred rate of size 
evolution (Mancini et al."2010'). 

There are several possible approaches to test the 
smallness of g: alaxie s at high redshifts. Recently, 
iSzomoru et al.l (|2010l ) used deep observations with the 
Wide Field Camera 3 (WFC3) instrument on board the 
Hubble Space Telescope (HST) to measure the size of a 
massive quiescent galaxy at 2: ~ 2 based on a new ap- 
proach (correcting the best-fit Sersic profile of the galaxy 
with the residual of the fit) to confirm the compactness 
of this massive galaxy at this redshift. 

The other method in order to check the effects of cos- 
mological redshift on the size/shape measurements is 
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to artificially transform nearby galaxies to higher red- 
shifts. Comparing derived parameters before and after 
redshifting provides a test for biases that may be in- 
troduced by degraded resolution and cosmological sur- 
face brightness dimming. This technique has been used 
in the past for different purposes, for inst ance assess- 
ing morphologies at higher redshifts (e.g., IPetty et al] 
2009"; 'Conselice et al.l 1201 ll: Fvan den Bergh et al.l [2QQ1 
Lisker et al. 2006; G iavalisco et al .1 Il996f ). Recently, 
Barden et al. (2008) used a set of ~ 100 local galaxies 
to study the cosmological redshifting effect on size and 
shape of galaxies at 0.1 < z < 1.1. They created new 
images from the best-fit single Sersic models of their in- 
put images, and then redshifted them to show that there 
is not any systematics on the size and morphological pa- 
rameters. However, nearby galaxies have signs of differ- 
ent sub-structures and low surface brightness features. 
Generating simulated galaxies with a comparable range 
of properties of galaxies and adding them into the blank 
sky background images is a practical test. However, these 
mock objects are simple cases compare to real objects 
and coul d be assumed as a lower limit on the systemat- 
ics (e.g.,|^jillo et al. 2006). 

It is also a common practice to measure the surface 
brightness profile of galaxies at high- 2; single compo- 
nent Sersic profile fitting. Therefore, it is assumed that 
for a consistent comparison of sizes at low and high- 
z, the profiles of nearby galaxies also should be mea- 
sured with the same method. However, as mentioned 
earlier, galaxies often consist of multiple components 
(e.g., bars, bulges, compact cores, spiral arms, etc,). 
In the local Universe, these sub-components are well- 
resolved and distinguishable in the photometric analysis 
of their structures. Therefore their surface brightness 
profiles may deviate from a single component model. 
It has been shown that using extra components in fit- 
ting surface brightness profiles of nearby galaxies bet- 
ter describes the underlying stellar distributions than 
using canonical single Sersic profile fitting (e.g., for El- 
liptical galaxies: [Fe rrarese et al. 1994; Laue r et "al]ll995l : 
iGraham et al.l 120031 : ^Huang e t al.l 1201 2. also see refer- 
ences therein) . Some authors have also shown that using 
single component Sersic profile fitting for nearby galaxies 
with more than one component might systen iatically bias 
sizes and morphological p arameters (e.g., iMeert et al.l 
[20TI iBernardi et al.ll20T2f ). 

Therefore, in this paper, we first investigate the bi- 
ases associated with estimating sizes of nearby galaxies 
using single Sersic profile fitting and its effect in their 
comparison to galaxies at high redshifts. These effects 
will also be tested against various types of galaxies (e.g., 
classifications according to their morphology, color, star- 
formation rate). We will explore the possible dependence 
of the systematics of sizes on the galaxies classifications 
and test alternate (two-component and nonparametric) 
methods. 

We also artificially redshift real images of nearby galax- 
ies (z ~ 0) to 2; = 1 in order to investigate the uncer- 
tainties of parameter measurements. We will use the 
resolution of HST WFC3 instrument, since images from 
this instrument are now being widel y used for studying 
galaxy structu r es at high redshifts (|Oesch et aTl l2010t 
Szomoru et a l.l 120121; IPatel et al.l 12012": 'Newman et al.l 
2012: ..Mosleh et al.l 120121 : Ivan de Sande et ah ^Oll , 



etc.). Moreover, for the sake of better statistics, we use 
a large sample of nearby galaxies (~ 1000 objects). 

Finally, we will use our robust size measurements to 
study the correlation of size and stellar masses of our 
nearby galaxies. Galaxies can be selected or classified by 
means of different methods or criteria, such as morphol- 
ogy, color and star- formation rate. We investigate the 
mass-size relation for different types of nearby galaxies at 
a wide range of stellar masses and test whether the selec- 
tion criteria could affect the mass-size relations. These 
relations provide a baseline for further studies at high 
redshifts. We will also compare the mass-size relations 
of galaxies after artificially redshifting them to 2: = 1 and 
examine if the robustness of galaxies mass-size relations 
at high redshifts. 

We explain our sample used in this study in section 
2. The size determination methods and their systematic 
offsets dit z ^ are explored in section 3. The stellar 
mass-size relations of nearby galaxies are studied in sec- 
tion 4. We describe the artificially redshifting procedure 
of galaxies to z = 1 and their sizes compare to z = 
in section 5. We discuss our results in section 6. The 
cosmological parameters adopted throughout this pa- 
per are Qm = 0.3, Qa = 0.7 and Hq = 70 km Mpc~^. 

2. DATA 

The sample of galaxies we use for this study is se- 
lected fro m the MP A-JHU SDSS DR7 (Kauffmann et_alj 
2003; Sali m et all [2007) which have sp ectroscopi c red- 
shifts for SDSS DR7 galaxies (jAbazaiian et al.l I2009D 
galaxies. The surface brightness limit for our sample is 
M50 ^ 23 mag arcsec^ with magnitude limit of r < 17.77. 
We initially select galaxies to have spectroscopic red- 
shifts within 0.01 < z < 0.02 and stellar masses of 
log (M^/M©) > 9. As we intend further to artificially 
redshift galaxies to 2: = 1, the imposed redshift limits 
are to avoid selecting galaxies where the SDSS PSF is 
broader than the WFC3 PSF at 2; = 1 (and provid ing 
sufficient sampling at high- 2;, see iBarden et al.lT2008L for 
more details), and also to avoid objects with very large 
apparent sizes. To reduce processing time, we further se- 
lect about 1000 galaxies randomly from this sub-sample 
(about one-third of galaxies in this mass and redshift 
range). We use SDSS r-band images for measuring their 
sizes at this low- 2;. 

We classify our sample into different sub-samples based 
on their color, morphology and specific star-formation 
rate (sSFR). The left panel of Figure [T] shows the distri- 
butions of all galaxies on color-magnitude diagram. The 
color and absolute magnitude are based on the New York 
University Value-A dded Galaxy Catalog (NYU-VAGC 
Blanton et al.|[2QQ5f ). Parallel to the red sequence distri- 
bution, we define the following line to separate galaxies 
into red and blue ones: 

{g-r)= 0.68 - 0.032(M^ + 20) (1) 

In order to classify galaxies based on their morphol ogy, 
we u sed the Galaxy Zoo catalog (GZl) (Lintott e t al.l 
&)Uh which is a morphological catalog of visually classi- 
fied SDSS galaxies. We classify galaxies into early-types 
and late- types based on debiasing fractio n of th e votes for 
each galaxy being dominant (see Lintott et aU (|2Qlll) for 
more details). We note that the classification are only 
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Figure 1. Left: The color-magnitude diagram of our sample. The solid line shows the separating cut to define red and blue galaxies. 
Middle: The same as left panel but galaxies are color-coded according to their morphological classifications, i.e. Early- type galaxies as 
red points and Late-type galaxies as blue points. The morphological classification is based on the Galaxy Zoo Catalog (GZl). Right: 
Distribution of galaxies specific star-formation rate (SSFR) as a function of their stellar mass; the solid line represents the separation cut 
at Log{SSFR) = -11. 



available for ~ 94% of our sample. The color magni- 
tude distributions of these early-types and late-types are 
shown in middle panel of Figure [J Early- type galaxies 
are indicated as red symbols and late- type ones as blue. 

Galaxies can also be selected by means of their specific 
star- formation rate (sSFR) (Brinchmann et al. 2004). In 
the right panel of Figure [U the distributions of sSFRs 
and stellar masses of galaxies are shown. We define the 
Log{sSFR) = —11 as a separating cut to split our sam- 
ple into star- forming and non star- forming ones. In sum- 
mary, we divide our galaxies by four criteria: (1) Mor- 
phology based on Galaxy Zoo visual galaxy classifica- 
tions (2) Color (3) Specific star formation rate (sSFR) 
(4) Sersic indices (based on smoothed profiles of galaxies 
at low- 2;, see Appendix B). 

We also need to take into account the effects of sample 
selections on the completeness. We follow S03 to apply 
volume corrections to our sample {Vmax method). We 
give each galaxy a weight which is proportional to the 
inverse of the maximum volume out to which it can be 
observed. However, as our sample is limited to redshift 
ranges of z = 0.01 — 0.02, all galaxies have equal weights 
and hence our sample is not biased by stellar mass 
incompleteness down to lO^M©. It is worth noting that 
the fraction of galaxies which were not morphologically 
classified by Galaxy Zoo is ~ 6% on average and hence 
the effects are negligible. 

3. SIZES AT Z = 

As mentioned earlier, the well-resolved profiles of some 
of nearby galaxies could exhibit non- Sersic structures. 
Consequently, this raises questions about the effects of 
these structures on measurements of nearby galaxy sizes 
(i.e., half-light radii) with single-component, analytical 
models. In the following, we employ several methods to 
measure sizes of our z = 0.01—0.02 galaxies. These meth- 
ods can be separated into two main categories: 'paramet- 
ric' i.e., measuring the half-light radius of galaxies us- 
ing best-fit two-dimensional analytical models and 'non- 
parametric' from their observed one-dimensional light 
profiles and measuring their total fiuxes as described be- 



low. 

3.1. Parametric Methods 

To quantify the structural properties of galaxies with 
parametric me thods, we use the GALFIT v3 model- 
ing software (Peng et al. 2010). GALFIT measures the 
shape and size of each galaxy by finding a best-fit para- 
metric model of its two-dimensional surface brightness 
profile. It generates a range of profile models which 
are convolved with the PSF of the galaxy image and 
determines the best-fit model by comparing models to 
galaxy light profile and minimizing the of the fit. 
GALFIT c an fit one or more analyti cal functions such 
as Sersic (iSersid 119631 : iSersid Il968f ). de Vaucouleurs 
(|de Vaucouleurslll948f ). etc., to a galaxy light profile. 

In the following we outline the procedure for using 
GALFIT and measuring galaxies structural parameters 
from Sersic models. We first created a postage stamp 
for each galaxy from SDSS {r-band) imaging frames 
(2048 X 1448 pixels and pixel scale of 0.396'')- The 
postage stamp should be large enough to contains enough 
background sky pixels. We initially set our postage 
stamps to have widths of at least 1800 pixels. However, 
as our galaxies have large apparent angular sizes and they 
might be located on different positions on SDSS frames, 
the postage stamp sizes varies a bit for each galaxy. Nev- 
ertheless, our defined box-size value creates a postage 
stamp for each galaxy > 10 times larger than appar- 
ent galaxy sizes. These are sufficient for leaving the sky 
background as a free parameter during fitting procedure. 

In order to detect and mask neighboring objects, 
we use SExtractor (jBertin fc ArnoutsI Il996f ). For 
SDSS r-band images, we use the following SExtrac- 
tor configuration parameters for detecting sources: 
DETECT.MINAREA = 10, DETECT.THRESH = 
1.5 and ANALYSIS.THRESH = 1.5 and 
DEBLEND_MINCONT = 0.095. In addition, 
we further smoothed out the mask map created by 
SExtractor to reduce plausible bias of sky background 
estimations by contribution from undetected low fiux 
regions around nearby sources. We also provide the 
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Figure 2. Simulation (I): The comparison between sizes of simulated galaxies (models with 'single' Sersic profiles) and their recovered 
sizes (using single Sersic component fittings in left panel and using Non-parametric method in right panel) and Sersic indices (middle panel) 
after adding them into empty regions of SDSS r-band images. As plots show, there is no systematics in recovery of parameters of single 
Sersic model galaxies for both methods. 
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Figure 3. Simulation (11): The comparison between sizes (i.e., half-light radii) of simulated galaxies (models with 'double' Sersic profiles) 
and their recovered sizes using single Sersic component fittings (right panel) and non-parametric method (left panel). This shows that sizes 
derived from single Sersic profile fitting are biased for true nearby two-component Sersic profile objects. We note that this simulation does 
not include noise, in order isolate the biases caused by intrinsically complex structures. 

initial parameters for GALFIT, such as half-light radius, 
magnitude, position angle and axis ratio derived from 
SExtractor and initially set the Sersic index to value of 
2. 

The SDSS photo pipeline generates a synthesized PSF 
image at the central position of each galaxy using a pub- 
lished tool of Read Atlas Images Q. We use this code and 
extracted PSF images in r-band SDSS for each galaxy, 
separately. The PSF images required by GALFIT to 
convolve model images during fitting procedure. 

3.1.1. Single Component Sersic Profiles 

Our first adopted parametric model for describing the 
galaxies surface brightness is the one component Sersic 
model. Single component Sersic profiles are widely used 

^ |http : //www . sdss . org/dr7/product s/iinages/read_psf . html | 



for determining galaxy structures and properties, espe- 
cially for high redshifts galaxies. The Sersic function 
describes the surface brightness of a galaxy at radius r 
as 

S(r) = See-^-[(^/^^)'^''-^] (2) 

where rg is the half-light radius and Hg is the surface 
brightness at Vq. The shape of the galaxy profile is de- 
termined by the Sersic index n, and the value of bn is 
coup led to n (see for more details in iGraham fc Driverl 
[2QQ5I) . 

3.1.2. Two- Component Sersic Profiles 

Although the single Sersic profile describes the sur- 
face brightness of gala xies over a large dynamic range 
remarkably well (e.g., iKormendv et al.l I2QQ 9). the de- 
parture from the simple models can be used for di- 
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agnosing the formation of the galaxies. Specificahy, 
nearby ehiptical galaxies tend to show either 'extra 
light' or 'missing-light' in their central regions, depend- 
ing on their luminosity and different empirical func- 
tions (e.g., "core-Sersic " or "Nuker" law) have been 
used and suggeste d to parame trize these distinct com- 
ponents (Ferrares e et al.l 119941; l^ uer et al. 2007, 1995; 
Graham et al... .2001 iCote et al.l t2006: .Hopkins et al., 
2009b[ ). However, as our sample consists of wide ranges 
of luminosities and morphologies, we use double Sersic 
profiles which allow a variety of possible inner and outer 
profiles for each object (see iTurner et al.l 120121 ). Our 
adopted multi-component model is described as: 

I](r) = Seie-''-[('-/'--)''"'-il +Se2e-''"^l('-/'-»^)''"'-il 

(3) 

To compute effective radii, we first analytically recon- 
structed the sum of the deconvolved circularized surface 
brightness profiles of two components from the best-fit 
parameters and then computed their total fluxes and con- 
sequently their half-light radii. 

3.2. N on- Parametric Method 

We test the results from these analytical models 
against an independent non-parametric method. The 
non-parametric technique do not rely on the previous 
assumptions about the structure of galaxies. It bene- 
fits from the galaxy observed curve of growth. In brief, 
the observed intensity profile of a typical galaxy is mea- 
sured through elliptical isophotal fitting and from that, 
the growth curve of galaxy fluxes is determined. This 
provides the radius at which the flux reaches to half of 
the total value. 

In details, in order to measure the half-light radius 
of galaxies from this method, we need to integrate the 
fluxes of galaxies at different radii and find the radius at 
which the flux reaches half their values. For this purpose, 
we first extract the observed surface brightness profile 
of galaxies using the IRAF task ELLIPSE (Jedrzeiewski 
[1987). This procedure measures fluxes in isophotal el- 
lipses over the galaxy image and therefore can make it 
possible to generate one-dimensional surface brightness 
profile of a galaxy. 

The accuracy of this method depends on the precise 
measurements of galaxies total fluxes. Therefore we mea- 
sure the fluxes out to ~ 400 arcsec from the galaxies 
center. However, the surface brightness of galaxies is low 
at the outer parts and hence it is very difficult to define 
the exact edges of galaxies. Therefore, for measuring the 
fluxes in the outer parts, we extrapolate the total light 
of galaxies beyond their petroR90 radius (i.e., a radius 
containing 90% of the petrosian flux derived from SDSS 
DR7). This is done by fitting one-dimensional Sersic pro- 
files to these outer regions. By integrating the light pro- 
files estimated from our best-fit models to infinity, the 
total fluxes in the outer regions are estimated. More- 
over, in this way, we also estimate the sky background 
for each galaxy as sky value is left as a free parameter 
during the fitting procedure. Then for each galaxy, we 
integrate fluxes at different radii up to radius smaller 
than petroR90 from the fluxes measured by ellipse fit- 
ting and add them to the fluxes estimated in the outer 
region. This gives the total fluxes of galaxies and hence 



leading to extract the radii at which contains half of the 
total fluxes, i.e., their non-parametric half-light radii. We 
note that for about less than 7% of the galaxies, the one 
dimensional fitting to the outer parts using petroR90 are 
not converged. Hence we used petroRbO which is radius 
contains 50% of the flux within petrosian magnitude. In 
order to check if the results depend on the choice of ra- 
dius for the rest of the sample, we repeated the procedure 
by fitting the outer parts of galaxies starting at smaller 
radii of i.e., petroRbO. The results were perfectly consis- 
tent for all galaxies. We note that we fixed the ellipticity 
{E) and the position angle (P-A.) of the ellipse isophotes 
to the values obtained from the best-fit of single Sersic 
parametric method. 

The sizes derived from non-parametric method are 
also needed to be corrected for point spread function 
(PSF) broadening and therefore we use the following 
relation according: R = A/(ri/2)^ — {rpsr)'^^ where ri/2 
and rpsF are the derived non-parametric half-light and 
PSF size, respectively. This correction is only important 
for very small galaxies (lesssiml kpc). It is also worth 
noting that that all sizes derived in this paper are 
circularized, using V^, in which a is the semimajor 
axis and b/a is th e axis ratio. This removes the effects 
of ellipticity fe.g.. ' Truiillo et al.l 120061 iFranx et al.ll2008l : 
>Williams et al 2010l ). 



3.3. Simulations (I) & (II) 

We perform simulations for testing our methods and 
procedures as follows. The first one is designed to test 
the reliability of the single component Sersic method 
and the non-parametric method for galaxies at z = 0. 
For that, we first generated single Sersic mock galaxies 
with random properties (magnitude, re, b/a) in similar 
range of values for real galaxies. We then add them into 
the empty regions of the r-band SDSS images and per- 
form our fits, using single Sersic profiles and the non- 
parametric method. The results are shown in Figure [2l 
In the left panel, the comparison between input and out- 
put sizes are shown using parametric method and the 
Sersic indices are compared in the middle panel. The 
output sizes derived by using the non-parametric method 
are also compared to their original sizes in the right panel 
of Figure [2l The sizes of galaxies can be recovered with- 
out any systematics with median differences less than 2% 
for both methods. There are also no systematic errors in 
the recovery of Sersic indices. This shows that for mock 
galaxies with single Sersic profiles, our procedure are ro- 
bust in recovery of their properties from SDSS images. 

As shown in Appendix A (Figure [TQ|). sizes of galaxies 
derived using single Sersic profile fitting can be biased, 
especially for massive early- type galaxies. This may be 
caused by the existence of additional component (s) or 
no n- Sersic light profiles. We have also shown that sizes 
derived with double Sersic component, are smaller that 
sizes from one component models. We test an idealized 
case using simulated two-component objects. For that, 
we first created a sample of 300 two-component Sersic 
galaxies such that each model galaxy has a central com- 
ponent with median half-light radius of ~ 1 kpc and 
the outer component with median size of ~ 3 kpc. We 
also assumed that the central components have larger 
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Figure 4. Top row: The stellar mass-size relation of late- type galaxies (left panel), and early- type galaxies (right panel). The individual 
galaxies are shown as small open gray circles and the blue and red filled circles are the median of the sizes in stellar mass bins. The 
solid blue and red lines are the best fits to the data. The mass-size relation from studies of Shen et al. (2003) and Guo et al. (2009) are 
also illustrated by dashed and dot-dashed lines, respectively. The best-fit relations are consistent with Shen et al. (2003), however, for 
early- type galaxies the relation fiattens below M <~ 4 x 10^^ Mq. The shaded gray regions show the physical sizes of PSFs in SDSS r-band 
images. Bottom row: The size dispersions as a function of stellar masses and their best-fits. The characteristic stellar masses, where the 
dispersions change significantly, are shown as vertical dotted lines. 



median Sersic indices than the outer-part components. 
For ah galaxies, the central components are ~ 0.6 mag- 
nitude fainter than the other one. These numbers are 
derived from the average results of the two-component 
fits to our real galaxies at z = 0.01 — 0.02. To ensure 
that we are testing only the effects of multi component 
galaxies, only the sky background levels are added to the 
images of these model galaxies without any additional 
noise or neighbouring objects. We then measure sizes of 
these two-component model objects using single Sersic 
profile fitting and the non-parametric method. The re- 
sults are shown in Figure [H As seen in the left panel, 
the sizes are recovered robustly with the non-parametric 
method. However, as shown in the right panel of Fig- 
ure [3l sizes from single Sersic fitting are biased (larger) 
compared to their input half-light radii, especially for 
large objects. This simplified test shows that sizes from 
single Sersic profile fitting can be biased for true two- 
component galaxies. Meert et al. (2012) use different as- 
sumptions for simulated SDSS galaxies and show the ex- 
istence of bias in recovered parameters by fitting a sin- 
gle Sersic profile to real two-component systems. Al- 
though our s ample are quite nearby (~ 45 — 85 Mpc), 
iBernardi et a l. (2012) show the same effect for the main 
SDSS sample at z ~ 0.1. Hence, using single Sersic sizes 
for local galaxies can introduce systematics in size anal- 
yses. 

Nevertheless, fitting correct models to nearby galax- 
ies is comphcated. Different authors use different mod- 
els to fit multi-component galaxies, e.g., traditional de- 



Vaucouleurs plus E xponential disk, Sersic +Exponential 
(Meert et al]|2012f ), double Sersic or eve n using multiple 
(3-4) Sersic profiles ([Huang et all 1201 2t ). It is also the 
case that not all of the galaxies (at wide ranges of stellar 
masses) need to be measured by multi-component mod- 
els (~ 77% are robustly fit with two-component models 
in this work). Therefore, for the rest of this study, we use 
our non-parametric sizes for these z ^ galaxies. Our 
simulations {I k II) demonstrate the robustness of our 
non-parametric method. In addition, due to the large 
angular sizes of our galaxies, the effects of the PSF on 
sizes from this method are negligible. 

It is worth nothing that fluxes used for estimating the 
stellar masses of SDSS galaxies are model dependent and 
hence these fluxes can be different from fluxes measured 
through the non-parametric method for each individual 
galaxy. Therefore, it is essential to correct the stellar 
masses according to the new flux measurements. We 
rescale the stellar mass of each galaxy by measuring the 
ratio between its non-parametric fiux and the one used 
for estimating its stellar mass from the MPA catalog. 

4. STELLAR MASS-SIZE RELATION AT Z = 

The stellar mass-size relation for the SDSS galaxies 
has been studied by SOS. They investigated this rela- 
tion for objects which are defined as early and late types 
according to their Sersic and concentration indices, and 
their relations have been widely used in literature. How- 
ever, it is argued that the half-light sizes used in SOS, 
which are from the NYU-VAGC catalog and based on 
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Figure 5. The stellar mass-size relation of galaxies classified by means of different criteria. In the right panel, blue, late- type (visually 
classified), star-forming galaxies, and low-Sersic index (n < 2.5) are compared and in the l eft panel re d, early- typ e (visually classified) , 
non star-forming, and n > 2.5 are compared. The stellar mass-size relation from studies of IShen et al .' (2003) and I Quo et al.l (j2009) are 
also illustrated by dashed and dot-dashed lines respectively. The points are the median size dispersions as a function of stellar mass and 
the lines represent the best-fits to these points. As this plot shows, the relation based on different methods of classification of galaxies are 
largely consistent, though blue galaxies lie above the other relations. 



Table 1 

The fitting results of the parameters in the Size-Mass relations. 



Sample a ^ j Mq ai (72 

Early-type -0.02 ib 0.07 1.25 ib 0.21 0.24 ± 0.73 10.67 ±0.20 0.74 ±0.07 -0.08 ± 0.24 

Red 0.04 ±0.05 0.80 ±0.12 -0.31 ± 0.47 10.53 ±0.13 0.75 ± 0.09 0.13 ±0.07 

Log{SSFR) < -II 0.01 ± 0.06 0.91 ±0.16 -0.05 ± 0.65 10.55 ±0.10 0.86 ±0.11 0.13 ±0.08 

Sersic > 2.5 0.09 ± 0.09 0.83 ± 0.21 -0.86 ± 0.90 10.53 ±0.16 0.75 ±0.15 0.14 ±0.09 

Late-type 0.06 ± 0.06 0.35 ±0.18 -0.20 ±0.55 10.60 ± 0.23 0.50 ± 0.04 0.16 ±0.10 

Blue 0.18 ±0.08 0.33 ±0.16 -1.40 ±0.75 10.32 ± 0.30 0.57 ± 0.05 0.05 ±0.12 

Log{SSFR)> -11 0.10 ±0.09 0.26 ±0.19 -0.74 ±0.83 10.20 ±0.21 0.66 ± 0.05 0.23 ±0.10 

Sersic < 2.5 0.12 ±0.08 0.27 ±0.16 -0.87 ±0.75 10.22 ±0.23 0.67 ± 0.05 0.24 ± 0.09 



Note. — The best-fit parameters for the stellar mass-size relation for different types of galaxies at z ~ (equation 4 and 5). 



one- dimensional (ID) sin gle Sersic fitting, could have 
been underestimated (e.g. lGuo et al.|[2QQ9l : lSimard et al.l 
I2Q11I ). We have also shown that using single Sersic fitting 
could bias the sizes of galaxies with high stellar masses. 
As the mass-size relation could depend on the fitting 
model employed (specifically at high-mass ends), our in- 
dependent non-parametric method for z = galaxies, 
should remove uncertainties due to model assumptions. 
Our sample consists of galaxies over a wide range of stel- 
lar mass (> lO^M©) and is suitable for investigating this 
relation. 

We first study the mass-size relation of our sub-samples 
based on morphological Galaxy Zoo classifications. The 
distribution of sizes versus stellar masses of late- type and 
early- type galaxies are illustrated in the top row of Figure 
m The median sizes in small bins of stellar masses for 
each sample are measured (blue and red circles), and it 
can be seen that sizes of both late- types and early- types 
show little or weak correlation with to mass up to ~ 
3— 4x lO^^M©; however the relation steepens beyond this 
stellar mass and is stronger for the early- type galaxies. 



For both types of galaxies, the relations seem to begin 
above specific stellar masses. 

To further quantify the correlations, we use the func- 
tional form employed for late-type galaxies in SOS (equa- 
tion 18) for both our late-type and early-type samples: 



7(M,/M0)"(l + M*/Mo 



,/3-a 



(4) 



Where, a, /3, 7 and Mq are free fitting parameters. 
This basically allows the relation to have two different 
slopes depending on the stellar mass range. The a and /3 
represent the slopes of the relation and the characteristic 
mass. Mo determine the stellar mass at which the slope 
of the relation changes. However, this relation is not very 
sensitive to the characteristic mass, Mq, therefore, this 
can be defined from the size dispersion relation as follows 
(equation 19 in SOS): 



O'lnR = 0-2 



((Ji - (72 ) 

+ (M*/Mo)2 



(5) 
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where ai and (72 are also free fitting parameters (repre- 
senting size dispersions at low and high masses) and Mq 
is the characteristic stellar mass at which cinR signifi- 
cantly changes. Size dispersions as a function of stellar 
mass for late-type and early-type galaxies are shown in 
bottom panels of Figure |4] (left and right panels, respec- 
tively). The best fits to the data points are shown as 
solid blue and red lines and the best-fit parameters are 
presented in Table 1. 

For late- type galaxies, the median size dispersions de- 
crease at stellar mass greater than ~ 4 x IO^^Mq con- 
sistent with 803. The mass-size relation for these galax- 
ies are also consistent with 803 (dashed line). The size 
dispersions for early-types also behave similarly and de- 
crease for massive galaxies above a characteristic mass 
around 4 x lO^^M©. However, due to low number of 
objects at these high mass bins, it is not clear how sig- 
nificant this effect is. 

The median sizes of early-type galaxies at a stellar 
mass range of log(M*/M0) ~ 10 — 11 are consistent 
with the 803 relation. However, at lower stellar masses 
(< 2 X lO^^M©), sizes are almost constant. Therefore, 
in this mass range, there is little or weak correlation be- 
tween stellar mass and size. For late- type galaxies the 
relation runs parallel at these masses but with larger 
sizes. The mass-size relations for galaxies with higher 
stellar masses (i.e., ^ 2 x lO^^M©) are steep for both late 
and early types. However, each sample exhibit different 
slopes and early-types have a steeper mass-size relation 
(see Table 1). 

We also present the mass-size relations of galaxies 
based on different sample definitions such as color, 8ersic 
indices and s8FR in Figure [5] in order to test the effects 
of these selections on the mass-size relation and defin- 
ing baselines for future studies based on different sam- 
ple classifications. Interestingly, the mass-size relations 
based on these classifications are consistent to the anal- 
ogous relations in Figure SI In the left panel of Figure [5l 
late-type galaxies are compared to the star-forming, blue 
and low 8ersic index galaxies. They are almost consis- 
tent, though the blue galaxies have larger sizes at stellar 
masses > lO^^M© compared to the others. This could 
be caused by excluding edge-on galaxies using the color 
criterion. We should note that the 8ersic indices are mea- 
sured from the degraded and smoothed 8D88 images of 
galaxies (8ee Appendix B, Figure [15]), hence removing 
biases from sub-structure. Nevertheless, it is interest- 
ing that the mass-size relation obtained for objects with 
8ersic indices of n < 2.5 are consistent with 803. The 
best fits to the mass-size relation are summarized in Ta- 
ble 1. 

The right panel of Figure [5] illustrates the compari- 
son between the mass-size relation of early-type galaxies 
and those with red colors, low s8FR and high 8ersic in- 
dices. The relations are also consistent with each other. 
For all samples, the relations are curved with a weak 
relation for galaxies below 4 x lO^^M©. The slopes of 
the mass-size relations at high mass ends (/3) for these 
red/quiescent /n > 2.5 galaxies are on average around 
~ 0.85, close to th e slope of early- type central galaxies 
in lGuo et al.l (|2009[ ). However, this slope is slightly larger 
for early- type galaxies. 

In general, we show that the stellar mass-size relations 
for both late and early type galaxies are curved with 



a steeper slope at higher stellar masses. The size 
dispersions bellow the characteristic masses are high 
but decrease above Mq. This is the case for all our 
studied samples. The stellar mass-size relations based 
on different definitions, such as color, s8FR, morphology 
are consistent with the scaling relations of late and early 
type galaxies. 

5. REDSHIFTING GALAXIES TO Z = 1 

In order to check whether cosmological effects and ob- 
servational uncertainties could affect size (and structural 
parameter) measurements of galaxies at high redshifts, 
we perform redshifting simulations of the low- 2; object. 
We use our sample of galaxies from 8D88 at 2: ~ to 
create artificially redshifted samples of galaxies resem- 
bling the same galaxies at z = 1 in H8T WFC3 images. 
Our redshi fting procedur e is similar to the method de- 
scribed by iBarden et aT] (^008) (FERENGI code) and 
we briefly describe it bellow. However, in order to take 
into account the effects of bandpass shifting, we only use 
8D88 r-band images as input and using WFC3- J125 
imag es from CANDEL8 DEEP DATA (Bouwens e t all 
120121 ) as output images inst ead of using t he k-correction 
method described in .Barden et al.l (j2008r ). WFC3 is the 
new near-IR instrument on board of HST and covering 
rest-frame optical wavelengths at z ~ 1 — 3. Hence it is 
suitable for this purpose. 

5.1. Method 

The first step in the redshifting procedure is to rebin 
the low-z images with pixel scale pi and redshift Zi to 
output images at redshift of (= 1 in this work) and 
pixel scale of Po by a factor of {3 as: 

J^o Po 

where D is the angular diameter distance, expressed 
as I) = -^Yqr^, and d is the luminosity distance. 

The next step is to apply cosmological surface bright- 
ness dimming at the rate of (1 + z)~^ in each rebinned 
pixel. By considering the fact that the absolute magni- 
tude of galaxies must be conserved, the total fluxes / of 
the input and output images must scale as: 



We note that it has been shown by sev eral studies (e.g., 
IBarden et all 120051 : iLabbe et al.l |2003") that the intrin- 
sic surface brightness of galaxies increases with redshift. 
Therefore, during our procedure of artificially redshifting 
our galaxies, we incorporate the surface brightness evo- 
lution and making the galaxies one magnitude brighter 
at z = 1 follow ing this equation M evo = xz -\- M and 
setting x = -l (IBarden et al.ll2008f ). 

It is important to reach the same resolution of the real 
data at high- 2;. Therefore, the next step is to correct 
the images to the appropriate P8F. This can be done by 
finding suitable kernels for convolving low- 2; images to 
reach the same P8F properties/shape at high- 2;. To do 
this, for each galaxy we require two P8Fs, i.e., its low- 2; 
and high-^ P8Fs. We use low-^ P8Fs from 8D88 (the 
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Figure 6. Example images of four spiral galaxies (left set of panels) and four elliptical galaxies (right set of panels). In each set, the 
left panels show the SDSS (r-band) postage stamp images (200^^ x 200^^) of galaxies at 0.01 < z < O.OlSand in the middle columns, we 
show their artificially redshifted (to z = 1) postage stamp (r^ x 7^^) images after adding to WFC3-Ji25 images and in the right columns 
represent the best fit single Sersic models of these redshifted galaxies. 

kernel depends on the input and output redshifts. 

After transforming images to the high- 2; resolution an 
pixel scale, the last step is to add background noise to 
the images. For this, we put galaxy images into ran- 
dom empty regions of the J125 band CANDELS DEEP 
images and then measure their structural parameters as 
described below. We note that, in order to check the 
effects of sky variations on galaxy property measure- 
ments, we repeated this step by inserting each redshifted 
galaxy into multiple empty regions. The final measured 
size/parameter for each object is the median of these 
seven realizations. 

The procedure to measure the structural properties of 
artificially redshifted gala xies (i.e., size a nd Se rsic index) 
is similar to that used in iMosleh et al.l (|2Q12f ) . In brief, 
we used GALFIT to find the best-fit single Sersic model 
for each galaxy. Neighbouring objects are detected by 
running SExtractor and masked during profile fitting. 
Initial parameter guesses, such as magnitude, half-light 
radius and axis ratio are provided from the SExtractor 
output. We used the median-stacked PSF from stars in 
the field. In Figure [H we show the SDSS postage stamp 
images of the galaxies (late-types in the left set of pan- 
els and early- types in the right set of panels) at low- 2; 
(the left columns) and after redshifting to 2; = 1 (mid- 
dle columns). The best-fit single Sersic models of these 
artificially redshifted galaxies at high- 2; are shown in the 
right columns. 

We perform two sets of simulations to test the size 
measurement accuracy in the J125 WFC3 images, and to 
check the procedure for artificially redshifting the galax- 
ies. These tests are described in Appendix B. We show 
that our redshifting method and size measurements at 
high- 2; are robust and can recover sizes and structural 
parameters of model galaxies without any systematics. 

However, as mentioned earlier, using single Sersic 
profile fitting for more complex galaxies in the nearby 
Universe potentially biases size estimates. This fact 
raises concerns about the sizes of galaxies at high 




—I I I L 



J I I I— 



—I I I I— 



0.0 0.5 1.0 
Log(re) (input) (kpc) 



1.5 



Figure 7. Sizes of simulated two-component 2; ~ galaxies which 
have been "redshifted" to 2; = 1 and remeasured with single- 
component Sersic profile fitting. The input and output sizes are 
consistent, indicating that sizes of multi-component galaxies can 
be reliably derived with single-component Sersic models at higher 
redshifts. 

ones we used for measuring sizes at z ~ 0) and for high- 2;, 
the median-stacked PSF which made from non-saturated 
stars in the J125 WFC3 images. Then, by transforma- 
tion of PSFs into Fourier space, finding their ratio and 
transforming the results back into spatial domain, we 
can find the convolution kernels required to reach the 
WFC3's J125 band PSF. Note that we calculate sepa- 
rately a transformation function for each galaxy as the 
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redshifts derived from single Sersic fitting. Therefore, it 
is also worth checking whether single Sersic profile fitting 
biases sizes of two-component objects at high redshifts. 
For that, we use the same simulated two-component 
model galaxies in section 3.3 (simulation (II), Figure [3]) 
and redshift them to z = 1. We measured their sizes 
after redshifting with single Sersic models. The results 
are shown in Figure [71 This shows that single Sersic 
profile fits of two-component galaxies at z = 1 provide 
reliable sizes, likely due to the smaller structures being 
washed out at high redshift. Therefore, traditional 
single Sersic surface brightness fitting robustly recovers 
sizes of our redshifted galaxies. 

5.2. Comparing to Sizes at z = 

In previous sections, we described how the sizes of our 
sample are measured reliably at both 2: ~ and z = 1. 
In this section, we compare sizes of galaxies before and 
after redshifting to z = 1. The comparison between low- 2; 
and high- 2; sizes for all galaxies are shown in left upper 
panel of Figure [8] and their median relative differences in 
small bins of sizes are shown in the left bottom panel. 
As can be seen, sizes before and after redshifting agree 
well and there are no systematics. 

There are also no biases if we split the sample into blue 
and red galaxies. Although the random scatter increases 
with size for red objects, there are no systematics, and on 
average sizes of these galaxies can be reliably recovered 
at high redshift. As mentioned in section 5.1, to check 
the effects of the sky background on properties of galax- 
ies at z = 1, we perform different realizations by insert- 
ing the galaxies into different random blank sky regions 
and remeasure their properties. The galaxy parameters 
at z = 1 are the median values of these repeated mea- 
surements and the error bars illustrate the la scatter. 
Galaxies are also color-coded according to their stellar 
masses. 

The results are the same using different galaxy classifi- 
cations. For instance, in Figure [HI the size comparison is 
shown for late-type (middle panel) and early-type galax- 
ies (right panel). Although the scatter increases for large 
and massive early- type galaxies, there is not any system- 
atic differences in their sizes. 

It is worth noting that for sizes at z = 0, we used non- 
parametric method and for galaxies at z = 1 single Sersic 
profile fitting is used. Using single Sersic profile fitting 
at z ~ brings systematic s when comparing; siz es before 
and after redshifting (e.g. IWeinzirl et all [201 1[ ). This is 
also the case for comparing Sersic indices, which tend 
to be overestimated at z ~ using single Sersic profile 
fitting. 

The fact that sizes of multi-component galaxies at 
z = 1 can be recovered robustly using single Sersic fitting 
can be explained by the resolution limit of images at high 
redshifts. The differences are mostly noticeable for mas- 
sive early- type galaxies. The bright center s of elliptical 
galaxies h ave typical sizes of < 1 kpc (e.g., Huang et al. 
I2QI2I : [Hop kins et al. 2009b. a) which is about the typical 
size of PSF FWHM of WFC3 images (- 1.2kpc) at z = 1. 
As a result, the inner components are smeared out and 
the galaxy profiles are dominated by the outer compo- 
nents. Therefore, using single Sersic fitting at this red- 
shift and resolution robustly recovers the true parame- 



ters (See Appendix B for additional tests which illustrate 
how degrading the resolution affects measured structural 
parameters of local galaxies). 

We have also checked whether the results after 
redshifting are sensitive to the signal-to- noise (S/N) of 
images. This has been tested by changing the S/N, 
either by adding noise to the SDSS r-band images before 
redshifting them or arbitrarily increasing the S /N of red- 
shifted objects. These tests did not show any systematic 
changes in the size of redshifted objects. Therefore, in 
general, sizes of galaxies at high redshift can be measured 
robustly using canonical single Sersic profile fitting as 
long as the physical resolution is not better than ~ 1 kpc. 

5.3. Stellar Mass-Size Relation after Redshifting 

It would now be interesting to examine how the 
stellar mass-size relations look like after redshifting to 
z = 1. In Figure [9l comparison of the mass-size relations 
before and after redshifting to 2; = 1 are illustrated. 
The relations for late-type galaxies are shown in the 
left panel, where the solid blue diamonds are the 
median sizes after redshifting in small mass bins and 
the dashed-three-dotted line is the best-fit to the data. 
The open circles are the median sizes at z = along 
with the solid black line as a best-fit (same as in Figure 
[4]). The size dispersions after redshifting also shown 
in the bottom left panel. The relations for early-type 
galaxies are shown in the right panel. As can be seen, 
the mass-size relations are consistent to their z = after 
redshifting. The poor constraints at the high stellar mass 
end are due to small number statistics; however, the 
results are consistent within uncertainties. This Figure 
shows that the stellar mass-size relations based on the 
single Sersic profile fitting at high redshifts are robust. 
Using different definitions for separating galaxies would 
result the same after redshifting. In Figure [9l we only 
present the relations for morphologically selected sample. 

6. DISCUSSION 

In this paper, we use a sample of about 1000 galax- 
ies at 0.01 < z < 0.02 from SDSS DR7 to study their 
sizes and stellar mass-size relations. We first investigate 
the robustness of size measurement methods for these 
nearby galaxies, using two main procedures for deter- 
mining sizes and structures: parametric methods (sin- 
gle and double component Sersic profile fittings) and a 
non- parametric method . In agreement with recent works 
(e.g. [Allen eFa l. 2006; Bernardi et al.[[2012[ : [Meert et al.[ 
[201 2! [Huang et al.ii2012[ ), the majority of galaxies in the 
nearby Universe are well fit with two component profiles. 
Comparing the sizes from the non-parametric method 
(and from double Sersic fits) to those from single Sersic 
fits shows the systematic overestimation of sizes from the 
single Sersic method. In particular, sizes and Sersic in- 
dices of early-type galaxies at the high stellar mass end 
tend to be overestimated using the single Sersic fitting 
approach. Non-Sersic profiles or substructures in nearby 
galaxies may to be the cause of this bias. We tested this 
by simulating two-component model galaxies and mea- 
suring their sizes using single Sersic fitting and reached 
the same conclusion. Using single Sersic profile fitting 
also overestimates the Sersic indices of these galaxies. 
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Figure 8. A comparison between sizes of galaxies at low redshift to their sizes measured after artificially redshifting to z = 1 for all 
galaxies (left panel), late- type galaxies (middle panel), and early- type galaxies (right panel). The errors are the standard deviation of their 
sizes measured at different positions (different realizations). Sizes of galaxies are recovered after redshifting without any systematics. 
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Figure 9. The stellar mass-size relations and size dispersions of late- type (left panels) and early- type (right panels) galaxies are compared 
before and after artificially redshifting to z = 1. The blue and red points are the the median sizes of galaxies in mass bins and the 
dashed-three dotted lines are their best fits. The mass-size relations are consistent with the relations at 2; = (open circles and solid lines). 
This further demonstrates that size-mass relations of galaxies are reliable at 2; = 1 using single Sersic profile fitting. 
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Therefore, we caution that relying on single Sersic fits 
can introduce biases for nearby, well-resolved galaxies. 

Stellar mass-size relations of z = galaxies from sur- 
veys like SDSS are often used as baselines for quanti- 
fying the evolution of higher-redshift galaxy sizes. Us- 
ing the nonparametric method and classifying our sam- 
ple thorough a number of frequently-employed criteria 
(size, color, morphology, and Sersic index), we have ex- 
plored the stellar mass-size relation of galaxies in the 
nearby universe down to a stellar mass of 10^ Mq. We 
show that the slope of the relation varies with mass 
for both late- type and early- type galaxies. The rela- 
tions flatten for gala xies below about 3 — 4 x 10 ^^M© 
(see also Figure 11 in lTurner et al.|[2Q12f ). Moreover, at 
these low stellar masses, the relations for both late-types 
and early-types run parallel but with smaller sizes for 
early- type objects. Above a characteristic stellar mass 
of ~ 3 — 4 X lO^^M©, the mass-size scaling relations get 
steeper with less scatter for both late- and early-types. 
However, the early- types have a significantly steeper re- 
lation than late- types. 

In S03, the mass-size relation for early- types is re- 
ported down to stellar a mass of ~ IO^^Mq. They in- 
dicated that faint ellipticals were missed in the analy- 
sis due to their type classifications based on concentra- 
tions and Sersic indices. However, they report tenta- 
tive evidence that the size-luminosity rel ation for faint 
red galaxies flattens at low masses. Grah am fc WorlevI 
(2008) showed that the size-luminosity relation of ellip- 
tical galaxies has a varying slope. Jan z fc Lisker (2008) 
also showed a different size luminosity relation for dwarf 
and giant early-type galaxies in the Virgo Cluster and il- 
lustrated that the relation has little to no dependence on 
luminosity on the faint end. Bernardi et aLj (j2012i) also 
pointed out the flattening of the early-type mass-size re- 
lation, in agreement with what we see in the right panel 
of Figure m With our \ow-z sample reaching lO^M©, the 
flattening of the relation at low masses is clearly seen. 

At the high stellar mass end (> 10^^ M©), we find that 
the sizes of the early-type galaxies tend to be slightly 
larger than the S03 relation. iGuo e t al. (2009) also found 
a similar trend for early-type central galaxies (dashed- 
dotted lines in Figure S]). We note that the early- types in 
our sample are morphologically selected and differ from 
the early- types in S03 (defined as n > 2.5). In addition, 
as mentioned earlier, sizes of high Sersic index galaxies 
maybe underestimated in NYU-VAGC. The number of 
our early-type galaxies in this stellar mass bin is low and 
provid es only weak constraints. However. . Bernardi et al.l 
(|2012l ) also show similar behaviour at this high-mass 
end. They note that the increase in the steepness of the 
mass-size relation for high-mass early-types could be due 
to BCGs (see also Bernardil 120091 : IBernardi et al.l l2007). 
They also pointed out that the steepness of the relation 
for early- types changes at these high stellar masses. 

However, it is still not clear how the massive early- type 
galaxies are connected to low mass ones (i.e.,< lO^^M©) 
and how the curvatu re of the mass-size relation arises 
for these galaxies. Gra ham fc W orlev (2008) argued that 
the curved size-luminosity relation for elliptical galaxies 
is expected from the assumption of varying profile shapes 
of these galaxies with luminosity a nd they are not dis- 
tinct types. However, Janz fc Lisken (|2008) find evidence 
for the different behaviour of faint and bright early-types 



(see also lToloba et al.ll2012f ). IBernardi etHI l\2012\ ) also 
pointed out that the curvature of the early-type scaling 
relations might arise from the presence of other com- 
ponents (e.g., a disk) with the bulges of these galaxies. 
On the other hand, the characteristic stellar masses men- 
tioned above are predicted by semi- ana lytical simulations 
for spheroids in IShankar et al.l ()2013f ) . They show that 
the physical processes behind the evolution of spheroid 
sizes are different below and above these masses, which 
might naturally explain the differing relations. 

For late- type galaxies, we find that the stellar mass- 
size relation is mostly consistent among our samples re- 
gardless of the exact definition. These relations are also 
consistent with the mass-size relation for late-types in 
S03. However, the mass-size relation for blue galaxies is 
somewhat offset to larger sizes and steeper than that de- 
rived for other "late-type" classifications. This is likely a 
co nsequence of t he str ong color-size relation pointed out 
bv iFranx et al.l (|2008f ). as well as the exclusion of red, 
edge-on spirals from the blue sample. The size-mass re- 
lationship among the "early-type" samples appears to be 
consistent regardless of the exact classification method 
used (elliptical, red, n > 2.5 and/or quiescent). 

Finally, we artificially redshifted our sample to re- 
semble z = 1 galaxies in WFC3 J125 band images and 
test the robustness of size and structural measurements 
at high redshifts. We remeasure sizes of galaxies with 
single Sersic profile fitting, a common method in the 
literature for high redshift galaxies. Our results show 
that using single Sersic profile fitting recovers the sizes 
of these redshifted galaxies without any systematics. 
Interestingly, this demonstrates that size measurements 
at high- 2; are robust, despite the single- Sersic models 
failing for nearby massive early- types. We further 
verified this with simulations, finding that once the 
small components of nearby two-component galaxies are 
smeared out at high z, single Sersic component fitting 
can adequately measure structural parameters. Image 
resolution is thus an important criterion for deciding 
whether to use single Sersic profiles. At a physical 
resolution < 1 kpc, where central bright components 
are well-resolved, overly simple models like single Sersic 
profiles can introduce biases, and a multi-component or 
non-parametric method should be used. 

7. SUMMARY 

We present the mass-size relation of a sample of nearby 
galaxies atz = 0.01 — 0.02, dividing the sample based on 
several common classifications. We examined different 
methods of size measurements in order to quantify the 
systematics associated in each method. We also artifi- 
cially redshifted these galaxies to z = 1 to test potential 
systematic effects on their size measurements at high red- 
shifts. From our results we find that: 

• Nearby early-type galaxies with masses > 2 x 
IO^^Mq are not well fit with single Sersic profiles. 
Two component fits and non-parametric methods 
appear to provide less-biased measurements. These 
methods give effective radii which are smaller than 
those measured with single Sersic fits. 

• The stellar mass-size relation of both late- type and 
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early- type galaxies are steep at high masses (~ 3 — 
4 X IO^^Mq) and flatten at low masses. 

• Although single-Sersic profile fits can be biased for 
nearby, well-resolved galaxies, they provide robust 
sizes at high redshifts. 

• The stellar mass-size relations of "spiral" and "el- 
liptical" galaxies are not particularly sensitive to 
the precise definition of these categories (color, 
Sersic index, morphology, sSFR), with the excep- 
tion of blue galaxies which follow a somewhat 
higher and steeper relation. 
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appendix 
appendix a 

FAILURE OF SINGLE SERSIC FITTING AT z = 

Surface brightness profiles of galaxies in the local Universe rarely conform to simple analytic models (e.g.. 'Allen e t al.l 
[2006; Simard et al. 2011). However, single Sersic profile fitting is widely used for measuring their structural parameters. 
In order to test whether or not using single Sersic fitting can bias the sizes of local galaxies with well-resolved profiles, 
we compare the half-light radii of our sample determined through different methods described in the text. In Figure 
[TQ| sizes of galaxies measured by using single Sersic profile fitting are compared to their sizes derived from the non- 
parametric method. In the left top panel, the comparison is shown for all galaxies and in the left bottom panel, the 
relative median differences of sizes as a function of single Sersic sizes are illustrated. As can be seen, for small galaxies, 
the median differences are small, however for large galaxies, the systematic differences reach to ~ 25%, i.e., sizes from 
single Sersic profile fitting are systematically larger than sizes from the non-parametric method for these galaxies. To 
diagnose the systematics, we show the comparison for blue and red galaxies separately in the middle and right panel 
of Figure [TOl respectively. This shows that for blue ones the systematics are less than ~ 10%, except for the large 
size end (~ 20%). However, for the red galaxies the systematic trend is significant and increases towards larger and 
massive objects (right bottom panel). 

We also compare half-light radii from one-component Sersic profile fitting and two-component Sersic profile fitting 
in Figure [TTJ As it shows, sizes from single Sersic fitting are on average larger for large and massive ones than sizes 
from two-component models. Specifically, for the red galaxies (right panels of Figure [TT]). there is a systematic bias 
toward larger sizes. 

The fact that sizes from single Sersic fitting are larger than sizes from the non-parametric method and two-component 
models, raises the question of how the half-light radii from one-component Sersic profile fitting could have been 
overestimated? As an example, a typical profile of an early-type galaxy is shown in Figure [121 In the left panel, the 
observed profile is shown as black circles and the non-parametric fit overplotted as a blue line. The red line is the 
best single Sersic fit to the galaxy. The single Sersic profile to the whole observed profile does not match completely. 
This can be seen from the extra-light in central regions of the residual profile which is illustrated in the lower left 
panel (green line). The half-light size from single Sersic fitting is illustrated by black diamond and is larger than 
the one derived from non-parametric method (black triangle). The sizes derived using the residual-corrected method 
(| Szomoru et al.l [20101 ) is also shown as a black star and as can be seen, this method also gives smaller size than the 
single Sersic profile fitting. 

The light profile of this galaxy can be described better by adopting two-components Sersic profiles. In the right 
upper panel of Figure [121 the two-component models and the total model are shown in dashed-dotted and solid red 
lines, respectively. The residual profile in the bottom right panel shows that this approach recovers most of the true 
profile of the galaxy. The half-light size derived from this method for this galaxy is consistent with the non-parametric 
size and hence smaller than single Sersic profile fitting. It is worth testing whether the choice of PSF could introduce 
uncertainties. For that, we remeasure sizes of this galaxy using a nearby non-saturated star as a PSF. This gives us 
the same results as before. Therefore, we conclude that the SDSS synthetic PSFs do not significantly affect. 

This basically shows that if massive galaxies are well-resolved or contain multiple components, structural measure- 
ments using a single analytical model could potentially be biased. We also used simulations to show this (see section 
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Figure 10. Top rows: The comparison between sizes of galaxies at 0.01 < z < 0.02 measured in two different ways, i.e., using single 
component Sersic profiles and non-parametric method, for all galaxies (left panel), blue galaxies (middle panel) and red ones (right panel). 
Bottom panels, show the relative differences between sizes of galaxies as a function of their one component Sersic sizes. The systematic 
differences between sizes of red galaxies increases up to about 40% and sizes based on single Sersic profiles are larger that non-parametric 
sizes. 
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Figure 11. The comparison between sizes of galaxies measured using one-component Sersic profile fitting and two-component Sersic surface 
brightness profiles (for ~ 77% of total sample). As it can be seen, galaxies which their profiles could be estimated by two-component Sersic 
profiles , have smaller two-component Sersic sizes compare to one-component Sersic sizes. 
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Figure 12. The observed profile (black open circles) of a typical early- type galaxy. In the left panel, the red line represents the best-fit 
one component Sersic profile and the blue line shows the best-fit one-dimensional Sersic fit to the outer part of the galaxy (> rpCtrogO for 
measuring non-parametric size). In the right panel, the solid red line represent the total best fit model from two-component models. The 
green lines in the bottom panels show the residuals from best-fits of one component Sersic profile. The size that derived using residual- 
corrected method (Szomoru et al. 2010) is shown by a filled star and non-parametric size is shown by filled triangle. The filled diamond 
shows the size of this galaxy using single component Sersic profile. This plot shows that profiles of galaxies at these very low redshift might 
be better explained by two components profiles. 

APPENDIX B 

HIGH REDSHIFT SIMULATIONS 

In order to check how weh we can recover galaxy properties in J125 WFC3 images, we perform simulations by gen- 
erating ~ 2300 synthetic simulated galaxies (assuming single Sersic surface brightness profile) with random properties 
within the following ranges: 20 < J125 < 26.5, 0.5 < n < 6.5 and < rg < 15 kpc^^i, convolving with the Ji25-band 
PSF. We further add sky background by inserting the simulated galaxy images into the empty regions of the real 
Ji25-band images, and then remeasuring their structural properties with the same procedure that we use for real 
galaxies. 

The results of the simulations are shown in Figure [131 Synthetic simulated galaxies are split into late- type (n < 2.5; 
left panel) and early-type (n > 2.5; right panel) according to their input Sersic indices. Then they are split into 
small bins over the size-magnitude plane. The relative differences between input and output sizes (i.e, A{re)/rein) in 
each small bin of size-magnitude distribution are measured and shaded accordingly. Then we overplot the distribution 
of artificially redshifted SDSS late-type and early-type galaxies on this size-magnitude plane, shown as blue and red 
points, respectively. Therefore, systematics in size measurements for each redshifted object can be estimated from this 
plot. 

As it can be seen in Figure [131 foi" late- type galaxies the systematic differences over the range of artificially redshifted 
SDSS galaxies are less than few percent. The systematic differences for early-type galaxies are also very small and 
only increase at the very faint magnitude end (i.e.Ji25 > 25). In general, comparing the distribution of artificially 
redshifted SDSS galaxies to the uncertainties in each bin shows that the systematics are expected to be very small 
(< 10% at most) for most of our sample. Therefore, we expect that our size measurement procedure at high- 2; recovers 
the properties of galaxies without introducing significant systematic biases. 

The next set of simulations is designed to check our redshifting procedure. For this purpose, we also create two- 
dimensional single Sersic model galaxies with a similar range of properties to the nearby SDSS galaxies. We assign 
them similar redshifts as our SDSS galaxy sample. Then we use our code to artificially redshift these mock galaxies 
to z = 1, insert them into J125 WFC3 images, and remeasure their properties using the method described in section 
5.1. The results are illustrated in Figure [Ml In the left panel, input sizes before redshifting and output sizes after 
redshifting are compared; in the right panel, the comparison of input and output Sersic indices is shown. The error 
bars come from the dispersion between different realizations (i.e., using different empty regions). This plot shows that 
the properties of these single Sersic model galaxies can be recovered after redshifting to 2; = 1 without any systematics 
and hence, our redshifting procedure works robustly. 

It is also discussed in the text that using single Sersic profile fitting likely measures the true structural parameters 
of galaxies at high redshifts. We verify this by simulating double-component galaxies and redshifting them to z = 1 
(simulation (II) & Figure [T]). In addition to these results, in order to test whether using low resolution images washes 
out the sub-components and changes the galaxies measured structural parameters, we remeasure galaxy sizes at 
z = using single Sersic profile fitting from their degraded images; i.e., images which are binned (by a factor of 
4) and Gaussian smoothed. Figure [151 shows the sizes (top row panels) and Sersic indices (bottom row panels) of 
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Figure 13. The relative differences between input and output sizes of synthetic simulated galaxies (simulation III) are measured in small 
bins over the (input) size-magnitude plane (left panel: n < 2.5 and right panel: n > 2.5). In each bin, the colors correspond to the median 
relative differences between recovered and input sizes of simulated galaxies. The red and blue points represent the artificially redshifted 
SDSS galaxies (Left panel: Late-type galaxies and Right panel: Early-type galaxies) on the size-magnitude plane. This shows that the 
systematics on size measurements of our galaxies are very small over their size-magnitude distributions. 




Figure 14. Artificial single Sersic model SDSS galaxies are created and then redshifted to z = 1 using our procedure. The left panel shows 
the comparison between their input sizes (at z ~ 0) to their sizes after artificially redshifted. The right panel illustrates the comparison 
between the input and output Sersic indices. The results indicate that our procedure recovers the parameters of Sersic model galaxies and 
there are no systematics on the measured sizes and Sersic indices after redshifting these mock galaxies. 

these galaxies after smearing. It can be seen that sizes and Sersic indices are smaher, especiahy for red galaxies after 
degrading. The results simply illustrate that the bright central parts of galaxies can bias the measurements of the 
structural properties of the galaxies in the nearby universe,when a single Sersic model is used. Also, it shows that 
resolution should be taken into account for structural parameter measurements. 
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